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In this work we study a class of leptophilic dark matter models, where the dark 
matter interacts with the standard model particles via the U{1) li-Lj gauge boson, to 
explain the excess in cosmic rays observed by ATIC and PAMELA experiments, 
and more recently by Fermi experiment. There are three types of U{1)l.-l. models: 
a) b) and c) Although ATIC or Fermi data is 

consistent with PAMELA data separately, ATIC and Fermi data do not agree with 
each other. We therefore aim to identify which of the three models can explain which 
data set better. We find that models a) and b) can give correct dark matter relic 
density and explain the ATIC and PAMELA data simultaneously recur to the Breit- 
Wigner enhancement. Whereas model c) with a larger Z' mass can explain Fermi 
and PAMELA data simultaneously. In all cases the model parameters are restricted 
to narrow regions. Future improved data will decide which set of data are correct 
and also help to decide the correct dark matter model. 

PACS numbers: 98.80. Cq ll.15.Tk 11.25.Hf 14.80.-j 



I. INTRODUCTION 



Recently, the ATIC and PPB-BETS balloon experiments have observed excess in the 
+ e~ energy spectrum between 300 and 800 GeV [l, 2]- The PAMELA collaboration has 
also reported excesses in the positron fraction from 10 to ~ 100 GeV, but shown no excess 
for the antiproton data j^, 4| compared with the prediction in the cosmic ray physics. These 
results are compatible with the previous HEAT and AMSOl experiments (e.g., p, la, LZJ) with 
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higher precision. Newly pubhshed result from Fermi also shows an excess at the e"*" + e~ 
energy spectrum above the background of the conventional cosmic ray model. However, it 

The excesses may be explained by astrophysical 



shows a softer spectrum than ATIC 



processes, for instance the nearby pulsars 



0,0, 



ll| or photon-cosmic ray interactions 12| . 



or due to annihilation or decay of dark matter (DM) partic 



hood 



predominately into leptons (e.g., [isl. 



16, 
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es in our Galactic neighbour- 



191]). Similar explanation 



can also account for the excess observed at Fermi by assuming a softer injection electron 
spectrum from these sources 20|. PAMELA indicate that the dark matter is hadrophobic 
or leptophilic [isl . In this work we show that gauged U{l)Li-Lj models proposed some time 
ago 21| in searching for simple Z' can naturally explain the excess in electron/positron spec- 
trum through the DM annihilation mechanism. Here Li and Lj are one of the three family 
lepton numbers with i ^ j. 

In constructing dark matter models, one should note that although ATIC or Fermi data 
is consistent with PAMELA data separately, ATIC and Fermi data do not agree with each 
other. Improved data are needed to decide which ones are correct and to distinguish different 
dark matter models. In the models we are considering, there are three different ways to gauge 
the lepton number differences: a) f^(l)Le-L^, t>) U{1)l^_l^, and c) U{1)l^_l^. Each of them 
has different features. Our aim in this work is to identify which of the three models can 
explain which data set better. 

The Z' in the gauged f/(l)i._i^ is leptophilic which only interacts with standard model 
(SM) leptons. If the Z' also interacts with DmQ [iqI ]. it will be the main mediator for 
DM annihilation with final product dominated by leptons offering possible explanation to 
electron/positron excess without anti-proton excess at PAMELA data. 

Another requirement for DM annihilation to account for the positron excess is that the 
annihilation rate {av) determining the electron/positron spectrum should be much larger 
than that determined from the usual thermal DM relic density. The enhancement factor, 
usually referred as the boost factor, is in the range of 100 ~ 1000. Several mechanisms 
have been proposed to produce a large boost factor, including_the DM substructures, DM 



nonthermal production mechanism 



the Sommerfeld effect 
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16 



23| and the Breit- 



Wigner resonance enhancement effect 2J]. The detailed calculation based on the N-body 
simulation shows that the boost factor from DM substructures is generally less than ~ 10 



25 



261]. In order for the Sommerfeld effect to be effective, the mediating particle needs 
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to be very light to allow long range interaction between DM particles. This mechanism 
in our case will result in a light Z' boson. However, there are very tight constraints on 
the coupling constant for such light Z' . Instead the Breit-Wigner resonance enhancement 
mechanism works very well in our model if the Z' mass is about two times of the dark matter 
mass. From the relic density and ATIC, Fermi and PAMELA data we can also constrain 
the f/(l)L^_i. charge of DM. There is no need to have very different [/(l)^-..^, charges of 



DM and SM leptons as in 18|. Further, the determined \J{V)i^._i^. gauge coupling strength 
is consistent with all constraints from LEP, {g — 2)^ and other data. 

The ATIC data also show a sharp falling at about 600 GeV in the electron/positron 
energy spectrum. We find such a feature needs substantial electron component as the DM 
annihilation products with the electron and positron pairs have a fixed energy from DM 
annihilation. Additional electron/positron energy spectrum with lower energy from sec- 
ondary decays will then help to enhance positron with lower energies. We find that the 
gauged U{X)Le-L^, (model a)) and \J{X)l^-l^ (model b)) can give excellent fit to the ATIC 
data while the f/(l)L^„L^ can not. On the contrary, Fermi shows a much softer electron 
spectrum, which does not favor the initial electron component. We find that the gauged 
f/(l)L^_L^ (model c)) gives an excellent fit to the Fermi, PAMELA and HESS data. 

The paper is organized as following: in Sec. II we give a brief introduction of the model. 
In Sec. HI we introduce the Breit-Wigner mechanism and the numerical results of relic 
density and boost factor in our model. Then we show the electron/positron spectrum of our 
model in Sec. IV. Finally we give discussions and conclusions in Sec. V. 

II. THE MODEL 



One of t 



anomalies 211] 



le following global symmetries in the SM can be gauged without gauge 



The gauge boson Z' resulting from one of the above models has the desired leptophilic 
couplings. At the tree-level the Z' only couples to one of the pairs e and e and r, and 
yU and r. We will use Y' to indicate the quantum numbers for one of the above three 
possibilities, Y' = Li — Lj. If the Z' in one of these models is the messenger mediating 
dark matter annihilation, the resulting final states are mainly leptonic states which can lead 



to electron/positron excess observed in cosmic rays. It is then desirable to have the Z' to 
couple to dark matter [l9|. We therefore introduce a new vector-like fermion with a 
non-trivial Y' number a. The reason for the dark matter being vector-like is to make sure 
that the theory does not have gauge anomaly required for self consistency. The Z' boson 
can develop a finite mass from spontaneous \J{\)i^^^i^. symmetry breaking of a scalar S with 
a non-trivial charge Y' = b. With the new particles Z', S and ip in the model, addition term 
Lnew have to be added to the Lagrangian beside the SM one Lsm with 

Lnew = -\z''^^'Z'^^ + Y,h''{-9%'Z'^)l + i'[Y{^^,-ag'Z'^)-m^]i; 

I 

+ iD^S)\D''S) + filS^S + XsiS^Sy + XsHiS^S)H^H , (1) 

where / is summed over the SM leptons. H is the usual SM Higgs doublet. 
The Z' coupling to fermions are given by 

L = -g'iaij^^ + la% - + ^i^'^Lu, - iyjYLiyj)Z'^ . (2) 

Note that the Z' coupling to leptons are flavor diagonal. There is no tree level flavor changing 
neutral current induced by Z'. 

After S and H develop non- vanishing vacuum expectation values (vev) vs and v, the 
physical components from S and H can be written as {vs + s)/V2 and {v + h)/V2, respec- 
tively. The non-zero vg will induce a non- zero Z' mass given by: m|, = h'^g''^v'g. A non-zero 
Vs together with the non-zero v will induce mixing between s and h with the mixing param- 
eter proportional to Xsrvvs- This mixing will change the masses and irth of s and h in 
the limit of without mixing. Since we will require m^/ to be much larger than the Z mass, 
this implies that is also much larger than m/^. The mixing will reduce the usual Higgs 
mass TTih. However, since the parameter Xsh is not fixed, if it is small enough the reduction 
in Higgs mass can be neglected. In any case, the effects of the mixing and also other term 
in the Higgs potential involving 5* and H will not affect our discussions in the following. We 
will not discuss them further here. 

The relic density of the dark matter is controlled by annihilation of ipil) Z'* lili + ViVi. 
The annihilation rate of dark matter av, with lepton masses neglected and summed over 
the two types of charged leptons and neutrinos, is given by 
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where v is the relative velocity of the two annihilating dark matter and s is the total dark 
matter pair energy squared in the center of mass frame. F^/ is the decay width of the Z' 
boson. If the Z' mass is below the ipil! threshold which we will assume, the dominant decay 
modes of Z' are Z' liU + z/jZ/j, and Vz' is given by, neglecting lepton masses 

r., = 4^,, . (4) 

In Eqs. (E]) and (jlj), we have assumed that there are only left-handed light neutrinos. If 
there are light right-handed neutrinos to pair up with left-handed neutrinos to form Dirac 
neutrinos, the factor 3 in these equations should be changed to 4. 



III. THE BREIT-WIGNER ENHANCEMENT AND BOOST FACTOR 

Since the relic density of DM is determined by the annihilation rate, the model param- 
eters are thus constrained. The same parameters will also determine the annihilation rate 
producing the electron/positron excess observed today, which requires a much larger annihi- 
lation rate. A boost factor in the range 100 to 1000 is necessary. We find that Breit-Wigner 
resonance enhancement mechanism works very well in our models if the Z' boson mass is 
about two times of the dark matter mass. 

The boost factor in this case comes from the fact that since the Z' mass mz' is close to 
two times of the dark matter mass m^, the annihilate rate is close to the resonant point and 
is very sensitive to the thermal kinetic energy of dark matter. To see this let us rewrite the 
annihilation rate into a pair of charged leptons as 

^ a'g'' 1 

where we have used the non-relativistic limit of s = 4m^ -|- m^f^, with 5 and 7 defined as 
m|, = 4m2 (1 - 5), and 7^ = r|,(l - 5)/ Ami. 

For thermal dark matter, the velocity f ^ is proportional to the thermal energy of dark 
matter. It is clear that for small enough 5 and 7, the annihilation rate is very sensitive to 
the thermal energy and therefore the thermal temperature T. At lower dark matter thermal 
energies, the annihilation rate is enhanced compared with that at higher temperature. This 
results in a very different picture of dark matter annihilation than the case for the usual 
non-resonant annihilation where the annihilation rate is not sensitive to dark matter thermal 
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energies. The annihilation process does not freeze out even after the usual "freeze out" time 
in the non-resonant annihilation case due to the enhanced annihilation rate at lower energies. 
To produce the observed dark matter relic density, the annihilation rate at zero temperature 
is required to be larger than the usual one, and therefor a boost factor. With appropriate S 
and 7, a large enough boost factor can be produced. 

For a detailed discussion, a precise form for the thermally averaged annihilation rate 



should be used which can be written as 

1 



24| 



{av) 



In 647r% 



n 



EQ 



Am?. 



a{s)\fsKi 



X 



ds. 



with 



n 



EQ 



27r2 X 



Ko(x) 



(6) 



(7) 



4m^ ■ av, 



where Qi is the internal degrees of freedom of DM particle which is equal to 4 for a vector 
fermion, Ki{x) and K2{x) are the modified Bessel functions of the second type. 




FIG. 1: The Breit-Wigner enhanced relative cross section (crv) / (av) x=20 as a function of time x. 



We plot in Fig. [T]the thermally averaged annihilation rate (av) as a function of cosmic 
time X = m^/T for several values of parameters 7 and 6. The annihilation rate when "freeze 
out" at X ~ 20 is adopted as a normalization to illustrate the enhancement effect today. 

To precisely determine the parameters we solve the standard Boltzmann equation of the 
decoupling process numerically. Fig. [2] shows an example of the evolution of DM abundance 
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FIG. 2: The evolution of DM abundance y as a function of x, compared with the equihbrium 
abundance Yeq oc x^^'^e~^ and the experimental measurements VLoh'^ = 0.1131 it 0.0034 [27I 1. 

Y = n/s with n and s the number density and entropy density respectively. The parameter 
ag''^ is adjusted to make sure that today's DM abundance is correct Y{x = 3 x 10^) = Yq = 
^0/^72.8 X 10S(m^/GeV) . 



IV. THE ELECTRON/POSITRON SPECTRUM 

We now discuss the DM annihilation produced electron/positron spectrum in the Ul^^l. 
DM models for all the three possibilities discussed before: a) — b) — L^; and c) 
L^ — Lr- In all these cases, 2/3 of the DM annihilate into charged lepton pairs and 1/3 into 
neutrino pairs, therefore the required boost factor should be 1.5 times of the case for DM 
only annihilate into charged lepton pairs. 

After the electron/positron pairs produced by DM annihilation they propagate diffusively 
in the Galaxy due to the scattering with random magnetic field 28|. The propagation 



processes in the Galaxy is calculated numerically in order to compare with data measured at 
the Earth. The interactions with interstellar medium, mainly the synchrotron radiation and 
inverse Compton scattering processes, will lead to energy losses of the primary electrons and 
positrons. In addition, the overall convection driven by the Galactic wind and reacceleration 
due to the interstellar shock will also affect the electron spectrum. In this work we solve the 



propagation equation numerically adopting the GALPROP package [29 1. 
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In Fig. [3] we show the model predictions on the /{e^ + e~) fraction and + e~ fluxes 
together with the observational data. The background is also calculated using GALPROP 
package 29] with the diffusion + convection model parameters developed in Ref. jl7|. Fol- 
lowing Ref. [30], to fit ATIC and PAMELA data we adopt the DM mass ~ 1 TeV and 
Merritt density profile. In addition for the U{1)l^~Lt (model c) we also plot the spectrum 
with DM mass 1.5 TeV to fit the Fermi result. A boost factor ~ 1200 (we have included the 
branching ratio into neutrinos which do not produce electron/positron excess)^, or equiva- 
lently {av) 3.6 x 10"^^ cm^ s~^, is found to give good description to the data 30|] (The 
cross section for l.STeV DM to fit Fermi requires the a slightly larger value of ~ 5.4 x 10"^^ 
cm'^s"^). Considering the errorbars of the ATIC data, we find that {av) = 2.7 — 4.5 x 10~^^ 
cm"^ s~^ (corresponding to a boost factor 900 — 1500) can be consistent with the observations. 
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FIG. 3: Ltft: positron fraction e'^ /{e^ + e ) predicted in the ?7(1)l-_l^. model compared with the 
observational data from PAMELA ^ and HEAT P, [sil] . Right: the total electron spectrum of the 



model, compared with observations of ATIC pj, PPB-BETS Fermi [8| and H.E.S.S. 



For our model scenarios, in the ) and b) one of the pairs is directly annihilated into 

electron/positron pair and another pair with secondary electron/positron pair. As can be 
seen from Fig. [3], these two cases predict a sharp falling of electron/positron in the spectrum 
at energy about 600 GeV and fit the ATIC data very well. We also note that the case b) is 
slightly favored over case a), but the difference is very small. 

For the case c), the electron/positron pairs come from secondary decay of muon pair and 
tauon pairs. If one normalizes the boost factor to PAMELA data, there is no sharp falling in 



^ Note that the boost factor in Ref. 30| is defined as (o'u)/3x 10^^ cni"^ s ^, instead of the one 
{(jv) / (av) x=20 as shown in Fig. [TJ 
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the electron/positron spectrum at around GOOGeV. This model cannot fit ATIC data well. 
Since this model has a softer spectrum, one wonders whether it can fit Fermi data. It is clear 
from Fig. [3] that for DM mass of ~ 1 TeV it can not fit the Fermi data either. However, we 
see from Fig. [3] that with a 1.5 TeV DM mass the soft electron spectrum can fit the Fermi 
data very well, and at the same time it can also account for the PAMELA data of positron 
ratio. 




FIG. 4: Constraints on the model parameters from the DM relic density and the cosmic ray 
data on the ^ — 5 plane. Dashed lines show the isolines with ag''^ adjusted to satisfy the DM 
relic density. Shaded region shows the allowed parameter ranges with additional PAMELA/ ATIC 
constraints, while the dot-dashed line shows the constraint from the best fitting cross section for 
PAMEL A/Fermi data. The a = 1 curve is given by requiring correct relic density. 

Finally combining the constraints on the model parameters from the relic density and 
PAMELA, ATIC/Fermi data we can derive the model parameters, which are shown in Fig. 
m By properly adjusting the ag'^, the whole parameter space of 7 — 5 can give the right 
relic density. However, to give the needed boost factor which can account for the cosmic 
ray electron/positron data, the parameters are further limited in a much narrower range, as 
shown by the shaded region in Fig. HJ which gives a fit to the PAMEL A/ ATIC data. For 
the model to fit PAMELA/Fermi data, we need a heavier DM mass and a bit larger cross 
section. The constraint is shown by the dot-dashed curve in Fig. HI We did not take the 
errors of Fermi data into account since they have very small statistic errors. We also plot 
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an a = 1 curve which gives correct rehc density while setting a = 1. The overlap between 
the a = 1 curve with the shaded region gives the correct relic density and boost factor to 
account for ATIC and PAMELA data with universal gauge couplings between dark matter 
and SM leptons. The allowed range of g' with a = 1 is 6.5 x 10^^ < g'^ < 9.2 x 10~^. Then 
we have 5.5 x IQ-^ < 7 < 7.5 x IQ-^ and 5 x IQ-^ < 5 < 1.2 x 10"^ for a = 1 to give 
correct relic density and boost factor. Similarly for the model to fit PAMELA/Fermi data 
we have g'"^ ^ 5.8 x 10~^, 7 ^ 4.8 x 10~^ and 6 ^ 3.5 x 10~^. In principle there should 
also be a narrow range satisfy the experimental constraints, which are not shown here. For 
a 7^ 1 case, one can extract a and g'"^ separately from Fig. [3] by the fact 7 ~ 3(y''^/12-7r in our 
models with a known ag'"^. An interesting thing to note is that the Y' charge of the dark 
matter does not need to be very different from the SM leptons. 



V. DISCUSSIONS AND CONCLUSIONS 



In this paper we studied a class of leptophilic dark matter model with U{l)L,~Lj gauge 
interactions to account for the recent cosmic ray results. We have shown that two of the 
anomaly free gauged [/(l)^-.^^ models, the gauged U{1)l^_l^ and f/(l)Le-L^ models, can 
naturally explain the excess in electron/positron spectrum at ATIC, while another model 
with U{l)Lf,-L^ gauge coupling can naturally account for the electron spectrum at Fermi. 
We have seen that in order to fit data the regions allowed for g'"^ and a are constrained 
to narrow regions. Plus the facts that the dark matter is required to be about 1 TeV to 
explain the sharp falling in electron/positron excess at ATIC or about 1.5 TeV to best fit 
the Fermi and HESS data, and Z' mass is required to be about two times of the dark matter 
mass to have a large boos factor via the Breit-Wigner enhancement mechanism, the models 
parameters are all determined in narrow regions. 

Since the model parameters are restricted to narrow regions, one has to check if the 
constraints on the models are compatible with other processes. A closely related process 
is direct dark matter search. The annihilation process, ip + ip + e~ , producing the 

relic dark matter density and the electron/positron excess, is related to the direct detection 
process by changing the interaction from s-channel to t-channel, i.e., + e —>■ tp + e. One 
may worry if the enhanced s-channel cross section, a large boost factor, will also lead to a 
large t-channel cross section resulting in conflict with direct detection results. This is not 
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the case. 

For ip + e ^ ip + e collision with mass of around 1 TeV ( or 1.5 TeV), the energy 
transferred to the electron is very small, at the eV order, far below the threshold energy of 
present detectors. However, if the dark matter particle ip kicks out a tightly bound electron, 
the energy transferred to the electron may be larger and detectable. But even this happens, 
the interaction cross section for the t-channel is still too small because in this case there is no 
Briet-Wigner enhancement, a{ip + e — ip + e) = {Qo? g''^ / 11)771*^ / m^, . Here m* indicates a free 
electron or a tightly bound electron. Using the bounds on the model parameters obtained 
before, we find that the cross section is of order (m*^/GeV^) x 10~^^cm^(2TeV/m2:')^; which 
is well below the sensitivities of the present detectors. 

Low energy processes such as g-2 of electron, muon or tauon, e — v collision data, and 
also high energy experiments at LEP have put stringent constraints on possible electron 
interaction with new gauge particles. However, in contributions from Z' are all proportional 
to g''^/m\, which is of order 10^^^/a(GeV~^)(2TeV/m^/)^ . If a is of order one there is no 
conflict with data. 

Direct production of Z' will be very difficult at LHC since at tree level the Z' only interacts 
directly with electron, muon/tauon. ILC will not be able to produce directly Z' either since 
its energy is below the Z' mass. But Z' can be produced at CLIC for models a) and b) 
where the e"'"e~ center of mass frame energy can be as high as 3 TeV. It may be interesting 
to scan the machine energy around the Z' resonant region to discover it. The decay widths 
of Z' in models considered here are of order 10 MeV. This requires a very careful scanning. 
For model c), a muon collider may provide a chance to produce the relevant Z' boson. 

Finally we point out that our study can be easily extended to a larger group of dark 
matter models that are leptophilic. Once we assume the Breit-Wigner enhancement takes 
effect to explain both the relic density and positron excess today the model parameters can 
be determined within a quite definite range, which should be similar to the present results. 

In summary, in the present work we have constructed simple dark matter models with 
gauged U{1)L^-Lj interactions. We have shown that the direct electron component in the 
annihilation final states are necessary to account for the sharp falling at the ATIC electron 
spectrum at ~ 600GeV. The U{1)l^^l^ and U{1)l^^l^ models with the dark matter mass 
around 1 TeV are in excellent agreement with ATIC data. However, to explain Fermi data 
which does not show the sharp falling indicated by ATIC data, one should avoid having 
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direct electron component. This fact favors ?7(l)L^_Lr niodel. All these models predict 
quite definite coupling constant, Z' mass and width, which are consistent with all present 
collider and other low energy data. 

At present ATIC or Fermi data is consistent with PAMELA data separately, but ATIC 
and Fermi data do not agree with each other. Although one can find models which fit one 
of the two data sets, i) PAMELA and ATIC, and ii) PAMELA and Fermi, it is obvious that 
these models cannot simultaneously all be correct. To finally decide which DM model is 
correct, experimental data have to give an unique set of data which can only be achieved by 
future improved experiments. Only by then one can decide the correct dark matter model. 
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